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Abstract 



(N 
(N 

, We consider a model for the injection-locked frequency divider, and study analytically the 

' locking onto rational multiples of the driving frequency. We provide explicit formulae for 

f-i I the width of the plateaux appearing in the devil's staircase structure of the lockings, and in 

particular show that the largest plateaux correspond to even integer values for the ratio of the 
frequency of the driving signal to the frequency of the output signal. Our results prove the 
experimental and numerical results available in the literature. 



^ '. 1 Introduction 

OO . In [27j . an electronic circuit known as the injection- locked frequency divider is studied experimen- 

\l ' tally, and the devil's staircase structure of the lockings is measured: when the ratio of the frequency 

uj of the driving signal to the frequency of the output signal is plotted versus uj, plateaux are 
' found for rational values of the ratio. In [5S], a model for the circuit is presented and numerically 

. investigated, and the results are shown to agree with the experiments. 

I In this paper, on the basis of the model introduced in [5S], we address the problem of explaining 

analytically the appearance of the plateaux of the devil's staircase. We aim to understand why 
the largest plateaux correspond to even integer values for the frequency ratio and, more generally, 
' how the widths of the plateaux depend on the particular values of the ratio. 

^ . From a qualitative point of view, the mechanism of locking can be illustrated as follows. For 

fixed driving frequency lo one considers the Poincare section at times t = 27rn/w, for integer n, and 
studies the dynamics on the attractor. This leads to a map which behaves as a diffeomorphism 
on the circle. Thus, based on the theory of such systems |2|, one expects that for lu/Hq close to a 
rational number one has locking. How close uj has to be to a rational multiple of fip depends on fj, 
and on the multiple itself: in the (oj, fi) parameter plane one has locking in wedge-shaped regions 
known as Arnold tongues. 

However, all the discussion above is purely qualitative. In particular, there remains the major 
problem of determining the map to which one should apply the theory. A quantitative constructive 
analysis is another matter, and requires taking into account the fine details of the equation and 
the explicit expression of the solution of the unperturbed equation: we carry this out in this 
paper. Our analysis is based on perturbation theory, which is implemented to all orders and 
proved to be convergent. This approach is particularly suited for quantitative estimates within 
any given accuracy (for which it has to be possible to go to arbitrarily high perturbation orders, 
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and to control the truncation errors). Furthermore, we think that a rigorous analysis ab initio, 
without introducing uncontrolled simplifications or approximations, can be of interest by itself. 
Indeed, although such simplifications can capture the essential features of the problem and allow 
a qualitative understanding of the physical phenomenon, it nonetheless remains unclear in general 
how far a simplified model can be expected to describe the original system faithfully. 

The conclusions of our analysis can be summarised as follows. The equation modelling the 
system can be viewed as a perturbation of order /i of a particular differential equation. In the 
absence of the perturbation, after a suitable change of variables, the system can be cast in the 
form of a Lienard equation x" + h{x) x' + k{x) — 0. Under suitable assumptions on h and fc, this 
admits a globally attracting limit cycle. Let ilo be the proper frequency of such a cycle, and let 
us denote by xa{t) ~ X{){D,Qt) the solution of the equation corresponding to the limit cycle, with 
the function being 27r-periodic in its argument. By also including the time direction, one can 
study the dynamics in the three-dimensional extended phase space {x,x' ,t), in which the limit 
cycle generates a topological cylinder. When the perturbation is switched on, the cylinder survives 
as an invariant manifold, slightly deformed with respect to the unperturbed case. This follows 
from general arguments related to the centre manifold theorem [5]. However the dynamics on the 
manifold strongly depends on the relation between the proper frequency SIq and the frequency 
(jj of the driving signal. If lo/Qq is irrational and satisfies some Diophantine condition (such as 
\uivi + rioJ^2| > 7(|j^i| + \v2\ + 1)^^ for all (t^i, V2) e and some positive constants 7, r), then one 
expects the output signal x{t) to be a quasi-periodic function with frequency vector uj = (a;, i7o), 
so that one has x{t) = X{u!t, Hot) = Xo(rio^) + 0(p), where X is a 27r-periodic function of both its 
arguments. In this case we say that the output frequency equals fio (of course, this is slightly 
improper terminology because Qq is only the frequency of the leading contribution to the output 
signal, and the latter is not even periodic). On the other hand, iiw/flo is close to a rational number 
p/q (resonance), then x{t) is periodic with frequency Q = pto/q (locking): hence the frequency Q 
of the output signal differs from Qq — even if it remains close to it — , because it is locked to the 
driving frequency uj. Thus, if one plots the ratio lu/Q versus lu one obtains the devil's staircase 
structure depicted in Figures 4 to 9 of ,27j . The locked solutions can be obtained analytically from 
the unperturbed periodic solutions by a mechanism similar to the subharmonic bifurcations that 
we have studied in previous papers [3l|5j. We stress, however, that, unlike the cases studied in the 
latter references, here, the unperturbed equation cannot be solved in closed form. This will yield 
extra technical difficulties, because we shall have to rely for our analysis on abstract symmetry 
properties of the solution, without the possibility of using explicit expressions. 



2 Model for the inject ion- locked frequency divider 

We consider the system of ordinary differential equations 

dVr dir 
C^=lL + f{Vc,t), L^ = -Rl^-Vc, (2.1) 

where L,C, R > are parameters, Vc and /l, the state variables, are the capacitor voltage and 
the inductor current, respectively, and 

f{Vc,t) = {A + Bsmnt)Vc{l-{Vc/VDDf), Vdd,A>0, B eR, 

is the (cubic approximation of the) driving point characteristic of the nonlinear resistor. The model 
p.l[) was introduced in [28j as a simplified description of the injection-locked frequency divider. 

By introducing the new variables u :— Vc/Vdd and v :— RIl/Vdd and rescaling time t — > 
Rt/L, (Pr]) becomes 

u' = av + <^(t)u{l-u'^) , v' = -u~v, (2.2) 



2 



where the prime denotes derivative with respect to time and we have set a — L/B?C, f3 — 
LA/RC, n = LB IRC, and = /? + ^sinwt, with cj = VlL/R. 

From now on, we shall consider the system ()2.2|) . with a > (3 > 1, and £ M. By setting 
a = u-\- V we obtain 

a' = [a-l)a + {^{t)-a)u-^(t)u^, u' ^ aa + {^{t) ~ a)u - ^(t)u^ , (2.3) 
which gives u" + [l - $(t) + 'i<^{t)u'^] u' + [{a - $(t)) u + $(t) w^] + <^'{t) {u^ - it) = 0, that is 
m" + (1 - /3 + 3/3u2) u' + [(a _ /3) u + /3 u^] + ^^f(u, u', i) = 0, (2.4) 

with 

*(m, u',t) = [u' (3u^ - l) sincji + {u^ - u) siwiot + uj {u^ ~ u) cosLot] . (2.5) 

For = 0, jMl) reduces to u" + (1 - /? + 3/3m2) + [(a - /3) u + /Jit^] = 0, which can be written 
as a Lienard equation 

u" + u' h{u) + k{u) ^Q, (2.6) 

with 

/i(m) = 1 - /3 + 3/3^t^ fc(u) = (a - /3) M + /3it3 = M (a - /3 + /3it2) . (2.7) 
For l|2.6p to have a unique limit cycle [51 [TH] , we require that 

l-/3<0, a-/3>0 =^ a>/3>l, 

which motivates our assumptions on the parameters a and (3. 

Consider the system described by the equation (|2.6p . with the functions h{u) and k{u) given 
by p.7p with a > (3 > 1. Such a system admits one and only one limit cycle encircling the origin 
|19| : cf. Figure [TJ Let Tq be the period of the solution uo(t) running on such a cycle. Denote by 
VLq = lix JTq the corresponding frequency: will be called the proper frequency of the system. 
Note that f2o depends only on the parameters a and [3. 

1.5 
1.0 
0.5 

0.0 

-0.5 
-1.0 



-0.5 0.0 0.5 

Figure 1: The limit cycle for a — 2.5, (3 = 2.0 and /i = 0. The proper frequency is Qo « 1.1434. 

The solution uo{t) is unique up to time translation. Fix the time origin so that u'q{0) — 0, 
Mo(0) > 0. Note that fixing the origin of time in such a way that Wq(0) — compels us to shift 
by some tg the time in the argument of the driving term in (|2.5p . i.e. '^(u, u' , t) must be replaced 
with \E'(m, u', t + to); cf. the analogous discussion in [TO] . 
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Lemma 1 The Fourier expansion of Uf){t) contains only the odd harmonics, i.e. 



{t) = e'"""' 



uo,, 



J2 e'^'^^'wo.. 



u odd 



Proof. The symmetry properties of p.6|) . more precisely the fact that h{~u) — h{u) and k{—u) 
—k{u), ensure that the periodic solution uo{t) satisfies the property 



uo{t + To/2) ^ -uo{t), 
and in turn this implies the result (compare the proof of Lemma 3.2 in [5 ) 



(2.8) 



Lemma 2 One has / dth{uQ{t)) > 0. 
Jq 

Proof. For a proof see [S] . 



Moreover the limit cycle is a global attractor [191 EI] > and it is uniformly hyperbolic [Ml IS] • 
Hence the cylinder it generates in the extended phase space persists, slightly deformed, as a global 
attractor for small perturbations [24l [8l EH [32] . This also means that the system described by 
the equation (|2.4p . at least for small values of /z, has one and only one attractor, and the latter 
attracts the whole phase space. However, the persistence of the attractor does not tell us whether 
the dynamics on the attractor is periodic or quasi-periodic; cf. [7| for an analogous discussion. In 
particular it does not imply that for Qq/u close to a resonance the dynamics remains periodic; cf. 
Figure H 




(a) 




(b) 



Figure 2: Examples of attractors for a = 2.5, P = 2.0 and fi — 0.1. For \uj — 4r2o| < 0.03 the motion 
is periodic: in (a), = 4Qo + 0.02. The black diamonds mark the four points where sinujt is zero and 
positive-going. For u — 4Qo + 0.2 the motion is quasi-periodic (b). Recall that flo ~ 1.1434. 

We note that for CIq/lo Diophantine the attractor is expected to become quasi-periodic, with 
the dynamics analytically conjugated to a Diophantine rotation with rotation vector {Qq,lu). In 
principle, this can be proved by KAM techniques [7], or with methods closer to those used in this 
paper [H HI [13 [M III] ■ 
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3 Framework for studying frequency locking 



Rescale time so that the driving term has period 27r, hence frequency 1, by setting r = cut. Then, 
by denoting with the dot the derivative with respect to rescaled time t, ()2.4|) gives 

M + - (l - /? + 3/3u^) u+ \(a- /3)u + Pu^] + /i^fu, m, r + To) = 0, (3.1) 

where tq = w^q, and we have defined ^{u,u,t) — [uj^^u{3u'^ — 1) sinr + lu~''^{u^ — u) sinr + 
u}~^{u^ — u) cost]. For /i = one has u + ^"-'^(1 — (3 + 3f3u^) ii + uj^^[{a — [3)u + j3 u^] ~ 0, which 
can be written as 

u+ -h{u)u+ ^k{u) = (3.2) 

UJ LJ^ 

which is of the form (|2.6p up to the rescaling of time. As an effect of the time rescaling, the 
frequency of the limit cycle for the system p.2p depends on as it is given by Via = Vla/uj. 

Remark 1 As the solution Uo{t) is analytic in t , the property 'ito(O) = means that we can write 
uq{t) = riT + O(t^), and hence Uo{t) = tq + riT^/2 + 0{t^), with — uq{0) and ri ~ iio(O). 

We want to show that if the frequency uj of the driving term is close to a rational multiple of 
the unperturbed proper frequency fJo of the system, that is w p^^j q for some p, g G IN relatively 
prime, then the frequency of the solution exactly equals qoj/p, that is a;/f2 = p/q- Such a 
phenomenon is known as frequency locking: the system is said to be locked into the resonance p:q. 
Let p = p/q e Q. For p, = 0, for any frequency uj of the driving term the proper frequency is 
= r^o — the system is decoupled from the perturbation — , so that if we fix w = pQo we obtain 
= Qq = Lj/ p. In terms of the rescaled variables, for which lj is replaced with tU = 1, the proper 
frequency becomes fig = l/p- For uj close to pilo write 

- = ^+e{p), (3.3) 

with e{iJ,) such that e(/Lt) — > as /i ^ 0. 

We look for periodic solutions for the full system (|2.4p . hence for solutions with period T = 
27Tp/uj (i.e the least common multiple of both 2tt/uj and 2T:p/ujq). In terms of the rescaled time t, 
the solution will have period 2ttp, hence frequency 1 /p. For = the system p.ip reduces to 

Ho{u, M, il) il + f{u) ii + g{u) = 0, (3.4) 

with 

f{u)^-^h{u), g{u) ^ ^-^k{u), 

which admits the periodic solution uo(t) such that uo(0) > 0, mo(0) = and mo(t + 2iTp) — uo(t). 
In other words the frequency of the limit cycle is 1/ p — q/p and the period is 2'Kp/q, i.e. uq{t) — 
U{t/p), with the function U being 27r-periodic. 
For /i 7^ we write 



e{p) ^eip + e2P + ... = 2_^ekP , (3.5) 

fe=i 

and, by inserting (|3.3p and (13. 5p into p.ip . we obtain the equation 

oo 

H{u, il, il, p) := Ho{u, it, il) + p^Hk{u, ii,T + tq) = 0, (3.6) 



fe=i 
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where 



Hi{u,u,t) = ei (1 - /3 + 3/5u2) it + -J- [(a - /3) M + /3 

P"0 



1 . o .^ . 1 . , ^ . 1 



+ — u (3m^ — l) sin T H — ^—^ (u^ — u) sin t H — — - (u^ — u) cos r, (3.7) 



H2{u,u,t) = £2(1-/3 + 3/3^2)?!+ f^+e?) [(a- /3)m + /3u3] 

+ eiii {iu^ — l) sin r H — ^ (u^ — m) sin r + ei (y? — u) cos r, (3.8) 

and so on. The shifting of time by tq = Lot^ in the driving term is due to the choice of the origin 
of time made according to Section [2l 

In the following sections we shall prove that for ^ small enough it is possible to choose e(^) as 
a function of Iq, in such a way that there exists a periodic solution of (|3.6p with period 27rp, i.e. 
with frequency When projected onto the {u,u) plane, such a solution is close enough to the 
unperturbed limit cycle (cf. for instance Figure [2|): the difference between them is of order /i. 

4 The linearised equation 

Write the unperturbed system (|3.4p as 

u = v, v — G{u,v), (4-1) 

with G{u,v) = -{pna)-Hl~ (3 + 3(3u^)v- ip^no)-^[{a -P)u + Pu^]. Let {uo{t),vo{t)) be the 
27rp-periodic solution of (|4.ip . which is uniquely determined by the conditions 

?io(0) = 0, uo(0) > 0. (4.2) 

The periodicity properties of uq{t) allow us to write 

Mo(r) - V e-^/^uo,., P=-, (4.3) 

^ q 



V odd 



as follows from Lemma [TJ Denote by 



^W2l(r) U)22(t) 

the Wronskian matrix of the system l|4.ip . that is the solution of the matrix equation 

where Gu and Gy denote derivatives with respect to u and v of G, and W2i{t) — wii(t), 'W22{t) — 

Wi2{t). 



Lemma 3 In \4-4^ one can set 



■wi2{t) := C2Uo{t), ■u;ii(t) ciuo(r) / dr' ^. ^ , (4.6) 



,ej^(-') 
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where F{t) is defined as 

F{t) -.^ [ dr' f{uo{r' j), f{u):^^{l-P + 3(3u')^^h{u), (4.7) 
Jo P^h P"o 

the constant f G (0,7rp) is chosen so that wii{0) ~ 0, and the constants C\ and ci are such that 
wii(O) = ^22(0) = 1. 

Proof. It can immediately be checked that {'Wi2{t),W22{t)), with wi2{t) defined as in (|4.6p and 
'u^22(''') = wi2{t), solves the linearised equation to l|4.ip . Then a second independent solution is of 
the form {wii{t),W2i{t)), with wii(t) given by (|4.6p and W2i{t) = wii{t); cf. [22], p. 122. In 
Appendix [Al we show that it is possible to choose f e (0, irp) in such a way that wii(O) = 0. The 
constants ci and C2 are chosen so that W{0) = 1. ■ 

Remark 2 With the notations of Remark{^ one has C2 = l/^i and Ci — — ri , so that C1C2 + 1 = 0. 
Note that uo(0) = 0, so that in (|4.6p the function u'ii(r) at r = is defined as the hmit 



,e-^(-') 



liniciMo(T) / dr' ..^ , (4.8) 



which is well defined; cf. Appendix [XI The same argument applies for r = 7rp, where uo{iTp) — 
uq{0) =0 — by (|2.8p . with the half-period To/2 becoming 7rp in terms of the rescaled variable. 

For any periodic function G we denote its average by (G) and set G = G — (G). Then /□ := 
if °uo) > (cf. Lemma [2|), so that we can write 



), F(t)= / dr' (/(zio(r')) - /o) , (4.9) 



where F{t), and hence e are well defined 27rp -periodic functions. 

Lemma 4 Given any periodic function P{t) and any real constant G 7^ there exists a periodic 
function Q{t), with the same period as P{t), and a constant D such that 

r dr' e^"'P(T') = D + cC-g(T). 
Jo 

One has D = -Q(0). 

Proof. Let P(t) be a periodic function of period T. Write 

P(r) = e""''-P., (4.10) 



where lo = 2ti/T. Then one has 

r dr' e^^'P(T') f e^"""'+^"' = ^ -P. 

so that, by setting 



„ G 4- iujv 



^^m^^Etvx— ^^=-E?^x— ' (4.11) 



G + iujv ^-^ G + illiv 



the assertion follows. 
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Lemma 5 There exist two 2tt p-periodic functions a{T) and 6(r) such that 

wiiir) = a(r) + e"^"^ &(r), w^ir) = ca(T), (4.12) 

for a suitable constant c. 

Proof. We cannot directly apply Lemma [3] because the function e~^^'^^ /uq{t) appearing in (|4.6p 
is singular. However we can proceed as follows. We write 



uHt) »)^Ou2(t)+77' 

so that the new integrand is smooth and it is given by e"-^""^ times a 27r p-periodic function. Hence, 
as long as < T < np, the integrand is bounded uniformly in 77, and we can apply Lebesgue's 
dominated convergence theorem, to write 

wii{t) = ciUo(t) lim / dr' .„ — . 

Then Lemma |4] gives 

wnir) = ciMr) Hm (e^"P(T, 77) - e^^P(f,77)) = ci lim {uo{r)e'^ P{r,v)) ~ Mr)e'' P{r), 

77— ►O r;^0 

where the function P{T,ri) is 27r/9 -periodic in r and P{f) — liiUjj^Q P(f,ri) is well defined. Note 
that —Uo{T)e'''^P{f) gives the function a(T) in (|4.12p . On the other hand, the function Wii(t) is 
also well defined, so that we can conclude that lim^^o {'Uo{T)e'^'^ P{t, r])) is well defined and smooth. 
As the function uo{t) P{t, rj) is periodic for any rj, the limit will also be periodic, and this defines 
the function 6(t) of (|4l^ . 

Comparing the expressions for wii{t) and Wi2(t) in (|4.6p . proportionality between the function 
wuir) and the periodic component of Wii{t) also follows. ■ 



Lemma 6 The Fourier expansions of the functions a(r) and &(t) in l[4.1S^ contain only the odd 
harmonics. 

Proof. Write ^0(1") according to (|4.3p . Then wi2(t) = ca(r) = cuo(t), so that the assertion 
follows trivially for a(r). Moreover the function e~^^'^'> /uq{t) involves even powers of functions 
containing only odd harmonics, so that it contains only even harmonics, and so does its integral as 
appearing in the definition (|4.6p of wii(r). Hence, by Lemma |4] and Lemma [SJ also &(t) in (|4.12p 
contains only the odd harmonics. ■ 

A straightforward calculation gives detM^(T) = — ciC2e~^*^'^-' = e~^^'^\ since C1C2 = —1, so that 

M/-i(r)=e^«f -""'fA. (4.13) 

We want to develop perturbation theory for a 27rp-periodic solution which continues the solution 
running on the unperturbed limit cycle when the perturbation is switched on. Therefore we write 

C30 

u{T)=UoiT)+Y,l^''MT), Ufe(T) = ^e^'^^/f?/fe,,, (4.14) 

k=l 1/6Z 

where uo{t) is the solution satisfying the conditions (|4.2p . Inserting (|4.14p into (|3.6p and expanding 
everything in powers of p, we obtain a sequence of recursive equations. In Sections [5] and [S] we 
shall consider in detail the first order. Higher order analysis and the issue of convergence will be 
discussed in Section [71 
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5 First order computations 

Let us also expand the initial conditions in fi: 

oo CX) 

L{Q)■.^u^uo{^) + ^^l''uk, u(0) := u = ^/wfc, (5.1) 



k=l k=l 



and set \E'i(r) = Hi{uq{t),vo{t),t + tq) — cf. (13.71) . We look for a solution {u{t),v{t)) which is 
analytic in /i, i.e. u{t) — uq{t) + iiui{t) + /i^M2(t) + . . . and u(t) = u(r) — cf. [31[51[Tnj for similar 
situations. Here we are interested in the dynamics on the attractor, hence in periodic solutions, 
but in principle we could also study the dynamics near the attractor, by looking for solutions of 
the form u{t) = [/(e^^^^, e'^^^, r), as in [I0l[l2l|4], with /i = /o + 0{y) and /s = 0(/z), and 
U{-, 27r/p-periodic in ip. To first order one has 



(5.2) 



and we can confine ourselves to the first component ui(r), since vi{t) — ui(t), 

ui{t) = wii{t)ui + wi2{t)vi + / dr'e^*^^ ^ [wi2(t)wii(t') - wn (r)wi2 (r')] ^'i(r). 

Jo 

which can be more conveniently written as 

m{T)^wii{T) (^Si-^dT'e^(^')u;i2(r')*i(T')) +u;i2(t) (^v, +^dr' e^(^')«;n(T')*i(r')) • 

The function e^^'^^wii(T)^'i(T) is periodic, while e-^'^'^'> wi2{t)'9 i{t) is given by e^""^ times a periodic 
function. Therefore we can write wii{t) and wi2(t) according to ()4.6I) . and set — cf. Lemma|4] — 

rdr'e^(^''a(T')*i(r') = e^«^Qi(T) - Qi(0), (5.3) 
Jo 

dT'e^("')e--^«"'&(T')*i(T') = tQo + Q2(r) - Q2(0), (5.4) 







for some periodic functions Qi(t) and Q2(t), and with 

Qo = {e^b'ifi). (5.5) 

Assume that we can choose the parameters in such a way that Qo ~ 0. Then we obtain 

ui{t) = a{T){ui+cQi{0)-cf'"^cQi{T)+cvi+e^°^cQi{T)-cQi{0) 

+ Q2(t) - Q2(0) ) + c-f"^biT) {ui + c Qi(0) - cf°^c Qi(t)) , (5.6) 

and if we want that (|5.6p describe a periodic function, the constant vi can assume any value, but 
we need 

ui = -cQi(0). (5.7) 

so that (|5.6p becomes 

ui(t) = a(r) (ci;i - cQi(O) + Q2(t) - Q2(0)) + c6(t)Qi(t), (5.8) 
where we have used that the function e-^"'^c Qi{t) appears twice but with opposite sign in (|5.6p . 



Remark 3 The constant vi is left undetermined, and we can fix it arbitrarily, say vi — 0, as we 
still have at our disposal the free parameter tq; cf. JT^, Section 2, for an analogous discussion. 

Therefore we can conclude that if Qo = then we can choose ui according to (|5.7p in such 
a way that up to first order there exists a periodic solution uq{t) + /iui(r) + 0{pi?). In the next 
Section we study in detail the condition Qg = 0. 
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6 Compatibility to first order 

Consider the equation Qo = 0, which can be written as 

exA + Bi{to) + B2{to) + B:,{to) = 0, 

where we have defined 



(6.1) 



A 



(l~(3 + 3f3ul) iiQ + — — (auo - (3ua + (3ul) 

P"0 



2Tip 



dTe^(")6(T) 



1 

— / dre^(")6(T) 
27rp Jo 



uo(t) h{uo{T)) + -7T-k{uo{T)) 

pih 

2 

uo{t) h{uo{T)) + —— /j(?io(r)) 
pilo 



(6.2) 



and 



Bi(ro) 
S2(ro) 
53(ro) 



dre^(^)6(T) 



/ dre-^(")6(r) 

27rp Jo 



1 

pilo 
1 



1 

27rp7o 



27rp 



dre^('''6(T) 



1 



7io(r) (3iig(r) - 1) 

2"o(t) (ug(T) - 1) 



pQo 



Mr) {ulir) - 1) 



sin(T + To), 
sin(T + To), 
cos(t + To). 



(6.3) 



Remark 4 Note that we can write {gj) as A ^ -L. J^'' ds Q{s) \^^h{U{s)) + ^k{U{s)) 

where U{s) and Q{s) are 2'K-periodic functions, with U{s) = uo{ps) and Q{s) = e^'''*-' b{ps). The 
function U{s) is the 2'K-periodic solution of the differential equation d^U /ds^ + ilQ^h{U) dU/ds + 
f2,^^fc(C/) = 0, and F{ps) — + ^ /o** ds' ft.([/(s')), so that the constant A is of the form 

A = A/p, with A independent of p. Hence if A ^ for some p £ Q then it is non-zero for all 
rational p =^ 0. 



By expanding sin(r + tq) = sinr costq + cosrsinro and cos(t + To) = cosrcosro — sinr sinro, 
we can rewrite (j6.3p as _Bj(ro) = Bn costq + Bi2 sinro for i = 1, 2, 3, where we have introduced the 
constants 



B21 
B22 



1 



— / dTe^(^)6(r) 



-i-7io(r) {iulir) - 1) 



2-Kp 



dre^'-^hir) 



1 

27rp Jq 

/ drc^(^)5(T) 

27rp Jo 



1 

pflo 
1 



27rp Jo 



27rp 



dTe-^(")5(T) 



1 



p^nl 



Mr) {iulir) - 1) 
2"o(t) {uI{t) - 1) 
Mr) {uUr) - 1) 



cos T, 



sinr = -pnoB32, 
cos T — pfloBsi. 



(6.4) 



By setting Di = - {Bn + B21 + B31) and D2 = - (B12 + -B22 + B32), (Ell) then becomes 

eiA — Sli(ro) := Di costq + D2 sinro- (6.5) 
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All constants Bij in (j6.4[) are given by the average of a suitable function which can be written as 
the product of a 27r/p -periodic function times a cosine or sine function. Consider explicitly the 
constant Bn; the other constants can be discussed in the same way. We write 

Bn = — drKMsiuT, with X(t) = V e'^^^/^if^ = V e^^'^^/'^ifa^, 

V even 

as follows from Lemmas |3] and El If we write sinr = X]ct=±i(^/^*)^*'^^' then 

i^ez, CT=±i 

2i/+(Tp=0 

The same argument applies to the other constants, so that we can conclude that the constants Bij 
can be different from zero only if p is an even integer. If we set p = p/ q this means q — 1 and 
p — 2n, n G IN. Hence for all rational p ^ 2E^ the first order compatibility equation ()6.5p gives 
eiA — 0, so that either A = and ei is arbitrary or A ^ and ei — 0. An explicit calculation 
(cf. Appendix [B]) shows that A ^ 0. Therefore for all resonances p:q, with p/q ^ 21M, frequency 
locking, if possible at all, can occur only for a range of frequencies of width at most /i^. 

The argument above does not imply that Di,D2 ^ for p/g £ 2W — in principle there could 
be cancellations in the sum (16. 6p . For any given resonance the non- vanishing of the constants 
Di and D2 can be checked numerically; for instance, when a = 5 and /3 = 4, for p/q = 2 one finds 
Di = 0.00735 and D2 = —0.04507 [6^. Therefore for p — 2n, n G frequency locking occurs for 
a range of frequencies of width of order fi around the value 2n. 



7 Higher order computations and convergence 

To extend the analysis of the previous sections to any perturbation order, we write the solution 
we are looking for as 



A:=0 



fc=0 



with (u(0),ti(0)) written according to (I5.1|l . Thus we find for all fc e M 



Uk{T) 

Vkir) 



W{t) 



\Vk 



+ / dr'T^-i(r') 





*fc(r') 



where 



Y /'i?fc'(u(r),u(T),r-t-To) 



.fc'=i 



(7.1) 



(7.2) 



(7.3) 



with Hh defined in (|3.6p . The notation [-Jfc for 5'fc(r) in (|7.3p means the following. In each 
term Hk', we expand u(t) and u{t) according to (|7.ip . and, by taking the Taylor series of the 
function Hk' , we keep all contributions proportional to /i*': we write the sum of these contributions 
as ^i^-^kir)- For instance one has ^'2(t) = H2{uo{t),vq{t),t + tq) + -i^Hi{uq{t),vq{t),t + 



ro)Mi(r) + ■^^Hi{uq{t),vq{t),t + to)ui{t), with ui(r) given by 

As in Section [5l we study only the equation for the first component, which is 



Uk( 



,{t) =wii{T)uu+wi2{T)vk+ / dr'e^^^') [u;i2(T)wii(r') - Wii(T)u;i2(T')]*fc(T'). (7.4) 
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The equation (|7.4[) for fc = 1 has been studied in Section[5l Here we want to show that the equation 
(|7.4p is well defined to any perturbation order k, and that it is possible to choose the constant Ek 
in (|3.5p so that it admits a periodic solution Uk (t) . 

The discussion proceeds as in Section [Sj once we note that each function Hk{u,u,T + tq) in 
(13. 6p contains a term (l — /3 + 3/3u^) -it + 2ek{p^o)^^ [(« — /?) w + Pu"^] , whereas all the other 
terms depend on the constants Sk', with k' strictly less than k. Therefore for fc S IN one has 

^kir)=ek{l-/3 + 3/3ul)uo{T) + ^[{a-l3)uoir)+f3ul{T)]+EkiT;To), (7.5) 

with the function Sfe(T;To) depending only on the constants Si, . . . ,£fe-i, besides the parameter tq 
and time r. Therefore to any perturbation order fc, in order to have a periodic solution, we need 

QkM-={c^b^'k) = 0, (7.6) 

and this can be obtained by requiring 

efeA = S)fe(ro), Dfe(ro) :=-(e^6Sfe(.;ro)), (7.7) 

with A defined as in ()6.2p . Since A ^ (as proved in Appendix IB|1 then wc can use (|7.7p to fix Ek 
as a function of tq. Defining the periodic functions Qk.iiT) and Qk,2iT) such that 

f AT'e^^-'^a{T')^k{r') = e^°^Q,,i(r) - Qfc,i(0), (7.8) 
Jo 



(7.9) 



/ dr'e^(^')e-/°^'6(T')*fc(T') = rQk.^ + QkA^) - QkA^), 
Jo 

choosing the constants Uk so that Uk + cQfe,i(0) = 0, and using (|7.6p . then (|7.4I) gives 

Mfe(r) = a(T) (c«fe - cQfe4(0) + QkA^) " Qfe,2(0)) + c6(T)Qfc,i(r), (7.10) 

with the constants Vk which will be fixed in the most convenient way (cf. Remark [J]). For instance 
we can set Vk — for all fc g M. 

We can make the perturbative analysis of the previous sections rigorous to all orders, by fol- 
lowing the strategy introduced in [16l [5], and hence study the convergence of the perturbation 
series. Alternatively, one could try to apply arguments based on the implicit function theorems. 
Typically, the latter would allow a simplification of the proof of existence of the periodic solutions, 
but would be less suitable for explicitly constructing the solutions themselves within any given 
accuracy; see the comments in [16]; therefore we follow the first method. Note that we are not 
confining ourselves to approximate analytical solutions, which could be unreliable because of the 
uncontrolled truncation of the series expansion. On the contrary we want also to settle the issue 
of convergence. In some sense this approach is complementary to that of [18' , where qualitative 
geometric methods are preferred to quantitative analytical ones. 

The study of the convergence of the series is standard, and it has been discussed extensively 
and in full detail in [16] for a similar situation. Thus, we only sketch how the argument proceeds. 

By expanding the functions u(r) and u(r) in Hk' (u(r), u(t), t + To) in (|7.3p according to (|7.ip . 
one sees that ^'^(t) can be expressed in terms of the functions Wfc'(r) with k' < k. On the other 
hand, by (|7.10p . the functions Ufc(T) are expressed in terms of the functions Qk,i{T) and Qk,2{T), 
which in turn are integrals of functions involving \['fc(r), and hence depend on Uk'^r) for k' < k. 

This means that we have recursive equations for the functions Uk (t) . By passing to the Fourier 
space, that is by expanding Ukir) = J2i^ex ^^"'^^^^k,!^, we obtain recursive equations for the Fourier 
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coefficients u^^u- We do not write them explicitly because the ensuing expressions are rather 
cumbersome, but one can easily work out the analytical expressions for the recursions by following 
the scheme that wG liRVG outlined-. Evcntucilly, we Cclh represent iik^u 

for fc > 1 and G Z, in terms 

of tree graphs, which can be studied with the techniques of |16j . 

We do not repeat the analysis here, but we instead just give the final result. To any order fc > 1 
one obtain the following bounds for the Fourier coefficients: \ukM\ < CiCj"^ and X^i/gz \uk,A < 
C3C'2~^i for suitable positive constants Ci,C2,C3, depending on p. This implies the convergence 
of the perturbation series l|7.ip for fi small enough, say for < Cj"^. 

8 Arnold tongues and devil's staircase 

We use the perturbative analysis, developed to all orders in the previous section, to study for which 
values of the driving frequency oj one has locking. We shall see that the analysis accounts for the 
devil's staircase structure found in ^28], for small values of the driving amplitude p,. 

Lemma 7 The functions Hk{u,u^T + tq) in Ii3.6\) are polynomials of odd order in {u,u) for all 
G IN. 

Proof. The function H{u, u, it, fi) given by (|3.H) is a polynomial of odd order in (u, u, ii). By writing 
H(u,u,u, n) as in (|3.6p . the only term containing il is the first one {k — 0), so that all the other 
terms are polynomials of odd order in {u, ii). ■ 

Lemma 8 For all k E ^ one has 

Mr) -EE e"'"^'e'"^"^^°^Uk,.,a (8.1) 

1/ odd |cr|<fc 

= E E e^'^^^e^'^'+'^^^fc.-.-^ (8-2) 

f odd \a\<k 

with the coefficients Uk^u. a and '^k,u,cr independent ofro. 

Proof. First of all note that if ^'fe(T) is of the form (|8.2p then Ufc(r) is also of the form (|8.H) . 
This can be proved as follows. For brevity, here and henceforth we say that Mfc(r) and ^'^(t) 
'contain only odd harmonics' if they are of the form (|8.ip and (|8.2p . respectively. The functions 
Qk,i{T) and Qk,2{T) are integrals of functions which are either periodic functions P(t) or of the 
form e-''"* times periodic functions P(t). In all cases the function P(t) is given by the product 
of three functions: two of these functions — one is either a(r) or &(r), the other one is ^'^(t) — 
contain odd harmonics, by Lemma[B]and by our assumption on ^'^(t), while the third one — e^^'''^ 
— contains only even harmonics. If we compare (|4.10p with (|4.11l) we see that the integral of a 
function e'~^'^ P{t) is of the form D + e'^'^Q{T), where Q{t) contains the same harmonics as P{t). 
Therefore both Qfc,i(T) and Qk,2{T) are periodic functions containing only even harmonics. Then, 
recall that Wfc(r) is given by (jT.lOp . We have already used the fact that the functions a{T) and 
&(t) contain only odd harmonics, so that we can conclude that, as claimed above, if \l/fe(r) is of 
the form (|8.2p then Uk is of the form (|8.ip . 

Then, the proof of the lemma proceeds by induction. Recall that for k = 1 one has 'I'i(t) = 
ifi(Mo(''"), ■uo(t), T + To), with Hi given by (|3.7p . so that, by Lemma [1] and Lemma [7l 4'i(r) is of 
the form (|8.ip . and, by the previous observation, the function mi(t) is also of the form (|8.2p . 



13 



By assuming that Ukir) is of the form (|8.ip for all k < k, then by Lemma [7] it also follows that 
^'j.(t), given by (18. 2p . is of the form ()8.2|) . Again by the observation at the beginning of the proof, 
it follows that Uj.(t) can be expressed as in ()8.2p . ■ 



Remark 5 If we expand Ukir) as a Fourier series, Ufe(T) = X^i^gz then i8.1\) implies 

v'e'£,\cr\<k 
qu +pcr— 1/ 

/n particular Uk^r) and ^'^(t) are polynomials of order k in Tq. 
Lemma 9 For a/Z fc e M one has 



Dfc(ro)=S)fc,o + — V V / dre-^/V^^^+^^^iffc,,.,, 
±:i tri Jo 



i/GZ (tGZ 
even 0<|cr|<fc 



for suitable tq- independent coefficients Kk^u.a, depending on £i, . . . ,ek-i, but not on Sk- 

Proof. The functions e^'"^-* and 6(r) in (|7.6p are periodic in r with period 2ttp = 2'Kp/q, and contain 
only even and odd harmonics, respectively, whereas ^'fe(T) is given by (|8.2p . By Lemma [U this 
yields that Qk^o = Qk,o{To) is of the form Qk,o = -ez E -ez /q dr e^''"/''e'-(-+^«)gfe.,,,, 

even \cr\<k 

for suitable coefficients Qk,v.<n which are independent of tq but depend on £i, . . . , e^. In particular 
the only contribution to Qkfl depending on Sk is of the form £kA — cf. (|7.7p — , so that we can 
write Qfc,o = SkA + Dfc,o(ei, ■ ■ ■ ,efe-i;To), for a suitable function Dfc^o(ei, ■ • • , £/c-i; To)- ■ 

By p.Sp and ()7.7p . and using Lemma [9l we can write 

£(m) -2)(ro,/i) ^X^/i'SfelTo), S)fc(ro) = J] e^'^^«I)fe,,, (8.3) 

fc=l o-GZ 

\a\<k 

for suitable coefficients Sfe.o-. For given w, for a periodic solution with period 2ttp to exist, we need 
that £(/i), defined according to (|3.3p . satisfy (|8.3p for some tq G [0, 27r). Therefore, by defining 

£max(/o) := max S)(ro,/i), emin(/3) := min I)(To,Ai), 

0<ro<27r 0<ro<27r 

and setting >V(p) = emax(p) — emin(p), such a periodic solution exists for all e(^) e VV(p). 

Lemma 10 Fix p — p/ q. One has 5)fc(ro) = Sfc^o for all k < q if p is even and for all k < 2q if p 
is odd. 

Proof. One can write EkA = !Dfc(To), with £)fc(ro) defined in LemmaO By comparing (|8.3p with 
the expression for £>fe(To) in Lemmali we see that Dk,a = 5^ E '^ez /g^''^ dr e"'^/Pe'''^''+'''>'> Kk,^^a 

^ V even 

for (T 7^ 0, so that one can have Sfc,^ ^ only if ap G 2W for some \a\ < k. Hence, if p — p/q with 
either even p and q > k or odd p and q > 2fc, one has 'I)k,a = 0. In other words, for fixed p = p/q 
one has ©^(to) = Sl^.o for all fc < g if p is even and for all fc < 2*7 if p is odd. ■ 



^k,cr, 



By Lemma \TU\ we can write in (|8.3p 



o-GZ 

0<|lT|<fc 



where the zero-mean function ©^(to) vanishes for /c < g if p is even and for k < 2q if p is odd. 
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Remark 6 The coefficient Dk,o does not contribute to W{p): when making the difference between 
£max o-nd Eniin Only Sfe(ro) plays a role. Therefore Lemma [TU\ implies that W(p) — 0{^i) only for 
p = 2n, n e ]N, >V(p) = ©(/i^) only for p = 2n - 1, n e K, W{p) = 0{p^) only for p = 2n/3, 
neK, Wip) = 0{p^) only for p ^ {2n ~ l)/2, n e = 0{p.^) only for p = 2n/5, neM, 

yV{p) = 0{p^) only for p ~ {2n — l)/3, n G ]N, and so on. In general, if p — p/ q with even p, then 
W{p) = Olp'i), while if p^p/q with odd p, then W{p) = 0(Ai^«). 

If we recall the definition (|3.3p of and we set 

l + pi2o£max(p) 1 + pS 2o £min (p) 

we obtain that for 

Wmin(p) < ^ < t^max(p) (8.5) 

there exists a periodic solution with period 2ttp (recall that p ~ p/q). In the (w, p) plane the region 
(|8.5p defines a distorted wedge with apex at uj — pQo on the real axis. 

Call Aw(p) = a;niax(p) — '^min(p) the range of frequencies around the value p^o, with p ^ p/q, 
for which there is frequency locking. Then 

ALj{2n/k) = 0(/), Auj{{2n + l)/k) = 0(p^'') (8.6) 

for all k,n € IN such that 2n/k and (2n + 1)/A:, respectively, are irreducible fractions. Indeed, 
Auj{p) is proportional to W(p), so that Dk,o does not contribute to the width of the plateau, but 
only to its 'centre'. In the {uj,p) plane the locking regions (Arnold tongues) 'emanate' from the 
values p^Iq, with p S Q. For p € 21N they are centred around the vertical passing through uj = pflo 
and for fixed p, have width 0{p). For all the other rational values of p, in general, they slightly 
bend away from the vertical: for fixed p the centre of the region is shifted of order p^ with respect 
to the value oj = p^lo/q, whereas the width is 0(/x') for even p and 0{p'^'^) for odd p. 



9 Conclusions and open problems 

The locking of oscillators onto subharmonics of the driving frequency (also called frequency de- 
multiplication) has been well known in electronics since the work of van der Pol and van der Mark 
[33] . In the (w, p) frequency-amplitude plane, the locking region occurs in distorted wedges (Arnold 
tongues) with apices corresponding to the rational values on the frequency axis. If one plots the 
ratio of the driver frequency to to the output frequency fl versus the driving frequency uj, one ob- 
tains a so-called devil's staircase, i.e. a self-similar fractal object, where the qualitative structure 
is replicated at a higher level of resolution, with plateaux corresponding to rational values of the 
ratio. 

The phase locking phenomenon, the existence of the Arnold tongues, and the devil's staircase 
picture have been proved rigorously in some mathematical models, such as the circle map [5], 
and studied numerically for several electronic circuits, such as the van der Pol equation 0^, the 
Josephson junction [lj[23l[30], the Chua circuit [29] among others. 

In this paper we have studied analytically the injection-locked frequency divider equation con- 
sidered in [28'. In particular we aimed to understand the devil's staircase picture, with the largest 
plateaux corresponding to integer resonances of even order, and to provide an algorithm to compute 
the width of the plateaux for small values of the driving amplitude p. 

The main result is summarised by (|8.6p . which gives the width of the Arnold tongues in terms 
of the driving amplitude p and of the resonances p : q. Note that the width of the tongues is 
narrower for resonances of higher order. 



15 



In most of the analytic discussions in the literature, one usually assumes that the unperturbed 
system is written in a very simple form — see for instance 'i8[. Of course, determining analytically 
the change of variables which puts the system into such a form can be very difficult in general, 
in principle as difficult as finding explicitly the solution itself. Hence, we have preferred to work 
directly with the original coordinates. Even if we have concentrated here on the injection- locked 
frequency divider equation, our analysis applies to any driven Lienard equation, of which the van 
der Pol equation is a particular type. The dynamics of the forced or driven van der Pol equation 
has been analytically investigated in [25l [26l El] • However, we could not rely on results existing in 
the literature, as we are interested in the exact structure of the Arnold tongues, which of course 
strongly depends on the particular form of the system under study. 

We have considered the model (|2.ip introduced in [28]. In particular we have taken the same 
driving term as in 28J, containing only one non-zero harmonic. In principle, one can consider more 
general functions, for instance any analytic periodic function, instead of the sine function. In that 
case the driving function contains all the harmonics; of course, by analyticity, the coefficients of 
the harmonics decay exponentially fast. Then one could ask how the analysis changes in such a 
case. From a technical point of view, there are no further complications. However, the conclusions 
about the devil's staircase structure are slightly different. For instance, the width of all plateaux 
becomes of order /i. This follows by the same arguments as given in Section [6l The analogues 
of the functions Bi{To) in (16. 3|) contain all the harmonics sin((T(T + tq)) and cos(cr(r + tq)), with 
0- G M, so that, when imposing the constraint 2iy + ap = Q in (j6.6p . one no longer has cr = ±1. 
On the contrary, one has cr G thus in general the constraint can be satisfied for all p G Q (by 
choosing ly appropriately), and so all the plateaux have width of order fi. However, the larger p 
and q in p = p/q are, the narrower the plateau is: indeed 2j/ + crp = (2i'q + (Tp)/q = requires 
2|z^|/|cr| — p/q^ hence, for very large values of p and q^ both v and a are very large, and hence 
the factors K21, contributing to Bn in (|6.6p are very small. This is consistent with the fact that 
the union of Arnold tongues form an open dense subset of the {i-u,p) plane, whose complement 
converges to full measure as /i — > [20,. So, an important observation is that large plateaux have 
not been found in [27\ for odd integer values because of the peculiar form of the driving term: they 
would appear by taking, for instance, a driving term involving also the harmonics with v = ±2. 

We have studied analytically the existence and properties of the periodic solution which contin- 
ues the unperturbed limit cycle when the perturbation is switched on. It would be interesting to 
prove analytically also that such a solution is attracting, for instance by determining the Lyapunov 
exponents or studying the more general solutions which move nearby and tend asymptotically to 
the attractor — for instance by following the strategy outlined in the first paragraph of Section [5l 

Another interesting problem to investigate analytically concerns the dynamics far away from 
the resonances, i.e. when the rotation vector (uj, f2o) satisfies some Diophantine condition such as 
the standard Diophantine condition mentioned in Section [1] — see also the comments in the last 
paragraph of Section [2] — or the weaker Bryuno condition [13l[T5]. Such values of w, in the devil's 
staircase picture, are complementary to those for which frequency locking occurs. 

The analysis we have performed is based on perturbation theory, and applies for p small enough. 
It would be interesting to investigate the locking diagram in the {i-u,p) plane for large values of p. 
It could be worthwhile to enquire further both analytically (for small values of p) and numerically 
(even for larger values of p) into the structure of the Arnold tongues in the {oj, p) plane. Work is 
underway concerning these problems |6j. 

Acknowledgements. We thank Giovanni Gallavotti for useful discussions, and Peter Kennedy for 
bringing this problem to our attention. We are also indebted to Henk Bruin and Freddy Dumortier 
for providing us with the references [S] and [34j . 
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A Well-posedness of the Wronskian matrix 



Let uq be the periodic solution of (|3.4p satisfying the conditions (14. 2p . Write uo{t) = ro + riT^/2 + 
rsrVS + 0(r4) — of. Rcniark[Tl 

Lemma 11 The function wii{t) in |^.^ is smooth. 

Proof. By deriving p.4p . one finds 

"0 + /(wo) "0 + /'(wo) + 5'("o) uo = 0, (A.l) 

where /' and g' are the derivatives of / and g with respect to their arguments, while the dots 
denote derivatives with respect to the time r. 

By computing (jA.ip at t = and using that 110(0) = 0, we find 

= uo{0) + /K(0)) MO) = 2r2 + /(ro) n. (A.2) 

In gJl) we can write F{t) = /g" dr' /(mo(0)) + ©(t^) = /(ro)T + ©(r^), so that e--^^^) = 1 - 
/(^■o)t + 0(t2). Q^jjgj. j^j^jjjj Qjjg j^g^g l/uo('^) = {rfr'^)^^ (l - 2r2T/ri + ©(r^)). Therefore 

the integrand in (j4.6p can be expanded as 

The term l/r^r^ produces a linear divergence, which is compensated by the function ?io(''') m front 
of the integral. The integral arising from the linear term inside the parentheses of (jA.3P would 
produce a logarithmic divergence (hence a divergence of the first derivative of wii(t)); however 
such a term is of the form — (2r2/ri + /(rg)) r = — rrj"^ (2r2 + /(t^o) ''i): which vanishes because 
of (jA.2p . Finally, the remaining part of the integrand arises from the terms of order in (|A.3p , 
and hence produces regular terms. This proves that the function ?i;ii(t) is smooth. ■ 



Lemma 12 There exists a unique r G {0,TTp) such that wii(O) = 0. 

Proof. One can write wii{t) in (|4.6p as wii{t) = ci'ito(''') (^(''") ~ ^(^))j where R{t) is a primitive 
of the function e-'^^^'> /uI{t), i.e. R{t) = r(T) := e-'^(^)/ug(T). 

The function r(r) is smooth and strictly positive for t G (0,7rp), and hence its primitive R{t) 
is strictly increasing for t e (0, 7rp). For all f e (0, 7rp) the function 



R{t, f) := J dr' r(r') = R{t) - R{f) 



(A.4) 



is smooth, and for all f G (0,7rp) one has lim^_>o+ i?(T, f) = —00 and lim^^^^p- _R(r, f) = +00, 
which imply that for all f 6 (0,7rp) the function i?(r, f) is strictly increasing in r from — cxd to 
+00. Now wii(t) = ci'Uo(t) (i?(T) — i?('f)) + ci e~-^^'^V'"o (''"), so that 

/ e--P'(^)\ 
?i;ii(0) = ci lim ilo{T)R{T) + - ciwo(0)i?(f). (A.5) 

T^o \ uo{t) J 

Lemma [Til shows that the limit in (IA.5P is well defined, so that we obtain wii(O) — provided 

1 / e-^^^'>\ 

^(^) = ^w7^ 1™ Mr) R{t) + . (A.6) 

Since R{f) is finite, by (|A.4p also the function R{t) is strictly increasing in r from — cxd to +00. 
Therefore (jA.Sp has one and only one solution f in (0, np). ■ 
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B Non- vanishing of the constant A 



Recall the definition (|6.2p of A. We can write e ^°'^h{T) — wii{t) — a{T) — wii(r) — ^uq{t), with 
7 = C2/C and uo[t) h{uo{T)) + 2(/90o)~^A:(uo(t)) = -pprioitolr) + uo(t) /i(uo(t))] — see (|4?6l) . 
(|4?T2l) and ([O]) — , so obtaining 

A---^/ ''dTe^(^)(u;ii(r)-7?io(r))[2pl7oMo(r)+Mo(T)/iK(T))]. (B.l) 



Lemma 13 One has 



2irp 



dTe^(^)uo(T) [2pnoUo{T) + uo(t) h{uo[T))] = 0. 



Proo/. By writing T{t) ^ e^(^), one has T{t) = f{uo{T))T{T) = h{uo{T)) T{t)/ p^o; cf. (l4Jl) . 
Hence J^iio [2pQouo + uo h{uo)] — p^o ( J'jt'^o + -^"o ) = P^o-^ (•^'"o) i so that 



27rp 



27rp 



dr e^(")uo(T) [2pnoMr) + Mr) Huoir))] = pfto J dr — (^(r)M2(T)) 

= pQa [Ti2Trp)ul{2Trp) - ^(0)^2(0)] - p^o {T{2Tip) - T{Q)) ul{0) = 0, 
where we have used that uq{t) is 27rp -periodic and uo(0) = 0. 

Because of Lemma [T31 (|B.1[) becomes 



(B.21 



Lemma 14 One has 



2Trp 



2-Kp 



dTe^(")wii(T) [2pnaUo{r) + Mr)h{uQ{T))] 



(B.3) 







dTe-'^^'^^u;ii(r) [2/9r2oWo(T) + uo(r) /i(mo(t))] = -27rpri/cir2o, 
with ri defined in Remark\^ 

Proof. By writing once more J-{t) = e^^'^\ we have 



dr 



where 



d 



QT 



dr 



(B.4) 



J- (wiiiio — wiiUoIt)) = — ^ {wiiW22 — W2iWi2{t)) = — !F det W = — J-e ^ = — = ri, 

C2 C2 C2 C2 

so that the integration of (|B.4|1 gives 

n2'JTp 

/ dTe^(-^^wn{T) [2pnoilo + M^) h{ua{T))] 

JO 

= pno [T{2TTp)uo{2TTp)wii{2Trp) - J^(0)wo(0)wii(0) + 27rpri] = 27rpripl7o, (B.5) 
where once more we have used that uo{2Trp) = uo{0) = 0. ■ 

By using Lemma [Ml in (jB.3|) we obtain A = — ripfio- Therefore A for any value p € Q. 
Note that the time rescaling implies that ri is of the form ri = (prio)^^fi, with fi independent of 
p, so that A = A/p, with A = — fi/f2o independent of p, consistently with Remark[4l 
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